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To understand the evolution of planetary systems, it is important to investigate planets in highly evolved stellar systems, 
and to explore the implications of their observed properties with respect to potential formation scenarios. Observations 
suggest the presence of giant planets in post-common-envelope binaries (PCEBs). A particularly well-studied system with 
planetary masses of 1.7 Mj and 7.0 Mj is NN Ser. We show here that a pure first-generation scenario where the planets 
form before the common envelope (CE) phase and the orbits evolve due to the changes in the gravitational potential is 
inconsistent with the current data. We propose a second-generation scenario where the planets are formed from the material 
that is ejected during the CE, which may naturally explain the observed planetary masses. In addition, hybrid scenarios 
where the planets form before the CE and evolve due to the accretion of the ejected gas appear as a realistic possibility. 
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1 Introduction 

While observational studies have predominantly focused 
on planetary systems around single stars, it is well-known 
that most of the stars are in binaries. It is therefore impor¬ 
tant to understand whether planets can form in binary sys¬ 
tems, which is suggested by an increasing amount of re¬ 
cent observations. The unevolved dG/dM binary Kepler 47 
harbors two planets with orbit al periods of 49.5 days and 
303.2 days ( Orosz et*^ 2012h . The binary s ystem Kepler 
16 ho sts the Saturn-sized planet Kepler 16b dPovle et al 


201 ih and Kepler 34/35 h osts a planet with 1/5 of Jupiter’s 


mass (IWelsh et al.L 120121) . From a theoretical point of view. 


To understand the evolution of planetary systems, it 
is particularly relevant to search for planets in highly 
evolved stellar systems, to determine if and under which 
conditions they survive different phases of stellar evolu¬ 
tion. The first exoplanets have been found around pul¬ 
sars, whic h provide accurate clocks for a precise measure¬ 
ment ( e.g.lKonacki & Wplszczan , 2003 ; Wolszczan & Frail , 
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formation of planets in gaseous disks formed by stellar 
winds in binary systems. 

An important class of highly evolved binaries where the 
presence of planets is considered are the post-common en¬ 
velope binaries (PCEBs). These are compact binary sys¬ 
tems of only a few solar radii consisting of a white 
dwarf and a low-mass main sequence star. For the PCEB 
system NN Ser, different series of mid-eclipse times 
have bee n obtained since the d i scove r y of the eclipses 


in 19 98 ( Beuemwmet^ , 2013 . 2010l : Brinkworth et al 


the for mation of planets in binary systems has been explored 
e.g. by lHaghighipour & Raymond (l2007h . 


20061: Haefner et all 2004). In 2010, the quality of the 
data was good enough to obtain a two-planet solution 
( Beuermann et al. . 2010h. which was shown to be dynam¬ 
ically stable jBeuermann et alil2013h . The most recent data 
suggest planetary masses of 1.7 Mj and 7.0 Mj with 
low eccent ricities and semi-major axis of ^ 3.3 AU and 
^ 5.4 AU (IBeuermann et al. . 2013b- The p lanetary solution 
has been confirmed by lMarsh et akl (1201 


1992t lYan et al.L 120131) . Recent obseryations haye re- 
yealed the presence of a p rotoplanetar y debri s disk around 
the p ulsar 4U 0142-1-61 ( Ertan et al. . l2007t IWang et dl 
2006b . Planet formation scenarios in h ighly e v olyed stel- 
lar systems haye been put forward by Perets ( 2010b and 
Tutukoy & Fedoroyal (12012b . considering for instance the 


At the same time, there is increasing eyidence against 
alternatiye explanations of the eclipsing time yariations. A 
frequently considered mechanism are periodic changes in 
the stellar structure related to magnetic actiyity, which could 
giye rise to regular changes in the star’s quadrupole mo¬ 
ment and there fore affect the gr ayitational coupling in the 
binary system (lApplegate[ll992b . Eor NN Ser, this mecha¬ 
nism has been excluded on energetic grounds, due to the low 
mass of the secondary star a nd the limited energy produc¬ 


tion (IBrinkworth et al.L l2006b . Indeed, one may expect that 


similar arguments can be made for many systems with low- 
mass companions (IZorotoyic & Schreiberil2013b . The pos¬ 
sibility of apsidal precession has also been excluded based 
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on recent data (IParsons et all l2014ll . leaving only the plan¬ 
etary hypothesis as a possible explanation. 

In this article, we will first review the main argu- 
me nts against a firs t -gener ation origin of the planets given 
by IVolschow et al.l (1201 4l) in section 2, and subsequently 
put forward a potential formation scenario for second- 
gener ation planets in section 3 (see Schleicher & Dreizleil 
20141) . Our main results and conclusions are summarized in 
section 4. 


2 The case against a first-generation origin 

To assess the feasibility of a first-generation scenario, we 
need to reconstruct the properties of the system before the 
common envelope phase. We will here discuss the case of 
initially spherical orbits, and refer to Volschow et al. (l2014l) 


for th e more general case. Following iBeuermann et al 


(120101) . we expect an initial size of the binary of about 
1.44 AU, corresponding to the typical size of a red gi¬ 
ant. The latter allows the system to enter the common en¬ 
velope phase during the later stages of the evolution. Us ¬ 
ing the binar y star e volution code by iHurlev et al. (2002), 
Mustill et al.l (l2013h have reconstructed the mass of the 
white dwarf progenitor to be between 1.875 Mq and 
2.25 Mq, assuming metallicities between 0.01 and 0.03. 
As the current mass of the white dwarf corresponds to 
Mwd = 0.535 Mq, we have a mass loss factor /i = 
■McurrentZ-^prev 0.3. The mass of the secondary corre¬ 
sponds to M 2 = 0.111 Mq, and plays only a minor role in 
the mass budget. In such a system, it is well-known that the 
planetary orbits are dynamically unstable, unless they ex¬ 
ceed at least three Hill radii. As a result, one expects semi¬ 


major axis of Tp 
we have 

GMini 


< = 


«/ = 


(3) 


2-/r-i’ 

with af the semi-major axis after the mass loss and r; the 
initial orbital radius. Assuming a typical expansion factor 


of 3, we therefore expect planetary orbits on scales of ^ 
10 AU, considerably larger than the observed ones. 

In these calculations, we have implicitly assumed that 
the mass loss occurs almost instantaneously. This is plau¬ 
sible, as most of the energy is released in the late stages 
of the common envelope phase, where the secondary has 
spiraled down to scales close to the central core, and the 
energy release occur s on very short timescales (see e.g. 


Kashi & Soken 1201 ih . However, the models for the com¬ 


mon envelope phase are still highly uncertain, and it is 
also conceivable that, at least for part of the evolution, 
the mass loss oc curs more gradually. Su ch a scenario has 
been explored by Portegies Zwar) ( 2013h for the planets in 
HU Aqua. In this case, the orbits of the planets will gradu¬ 
ally adjust according to the decreasing mass of the system. 

From the conservation of angular momentum, one can 
then show that rf = with ry and the hnal and 

initial radius of the spherical orbit. With the expected mass 
loss factor p, ^ 0.3, we therefore hnd an orbital increase 
by a factor of ^ 3, suggesting that the observed planetary 
orbits should be 10 AU, while in fact the planets have 
been detected with semi-major axes of 3.3 AU and 5.4 AU. 

Additional arguments can be made based on the ec- 
centricities of the pla nets. In the orbital solution of 


Beuermann et akl (120131) . these are given as ei = 0.144 and 


€2 = 0.222 for the outer and inner planet. In the case of 
an instantaneous mass loss, o ne would expect a sign ificant 
increase of the eccentricity as dVolschow et an. l20l4 
e/ = p-i-l. (4) 

For mass loss factors of p ^ 0.5, we therefore expect ec- 
centricities o f order 1. Similar results have been reported by 


Mustill et al.1 (120131) based on a dynamical analysis. 


2.5 AU. Assuming Keplerian rotation. 


( 1 ) 


where Mini is the initial mass of the system. Due to the 
mass loss of the system, the gravitational binding energy 
will decrease, so the orbits of the planets will widen or even 
become unbound. Considering that the escape velocity is 
given as 

2 _ 2GM 

^esc ; G) 


with M = fiMini the central mass after the ejection, the 
planets will become unbound for a mass loss factor p < 0.5. 
As mentioned above, a mass loss of p 0.3 is expected 
for NN Ser, and it seems likely that the previous planets 
became gravitationally unbound. Even for more moderate 
mass loss es, one can show that t he planetary orbits would 
evolve as (IVolschow et al.Ll2014l) 


To avoid these constraints, gas drag forces are neces¬ 
sary to decrease the planetary velocities by friction effects 
to keep them gravitati onally bound and at low eccentric¬ 
ities. For this purpose, Volschow et alJ ( 20141) considered 
both a spherically symmetric outflow assuming an exponen¬ 
tial density prohle, as well as a disk-shaped ejection mecha¬ 
nism in the orbital plane of the binary system. In the spher¬ 
ically symmetric case, their calculation shows that the drag 
forces is clearly subdominant, with a typical ratio of 10“^° 
compared to the gravitational force. For a disk-shaped out¬ 
flow, there can be some initial drag inwards, which is com¬ 
pensated during the subsequent evolution, where the drag 
forces are reduced during the expansion of the disk, and 
even an outward acceleration may occur when the disk ma¬ 
terial moves beyond the planetary orbit. In fact, in most of 
the simulations the planets are ejected, and it appears un¬ 
likely that the observed low-eccentricity orbits could be pro¬ 
duced. A potential way to avoid the problem could be an 
angular momentum exchange with an inhomogeneous gas 
distribution, which could however easily produce an infall 
of the planets onto the star. Alternatively, one could increase 
the effect of the gas drag and the resulting friction by main¬ 
taining high gas densities near the planets for several orbital 
periods. The latter can be the case if some material is ejected 
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but remains gravitationally bound, leading to the formation 
of a fail-back disk. In such a case, the gas will be accreted 
onto the planets and affect their dynamical evolution. In the 
next section, we will consider the formation of such a fall¬ 
back disk and the formation of new planets via gravitational 
instabilities. These planets may either form from seeds cre¬ 
ated by the gravitational instability, or potentially also from 
the previous planets, which may become the cores of the 
new ones. In the first case, we expect a truely second genera¬ 
tion, while the second case corresponds to a hybrid scenario 
between first- and second-generation planet formation. 


at a radius of 0.9 R©. In general, one may therefore ex¬ 
pect the mass in the envelope to be ejected, but one needs 
to determine which fraction of the mass becomes truely un¬ 
bound. The latter requires to also determine the velocities 
of the ejected material. Clearly, the latter can be pursued 
only in an app roximate fashion with in an analytical frame¬ 
work (see e.g. lKashi & Sokea 1201 ih . while numerical sim¬ 
ulations may not be able to add ress the later stages of the 
evolution (cf. lPassv et al.Ll2012h . 


3 Planet formation from the ejecta of 
common envelopes 

In this section, we describe the formation of a fail-back disk 
after the common-envelope phase, as well as the formation 
of giant planets via gravitatio nal instabilities as outlined by 
Schleicher & Dreizler ( 2014 ). 


3.1 Ejection and the formation of a fail-back disk 


For the ejection event, w e adopt the model oflKashi & Soker 
(|201 ih as employed by Schleicher & Dreizled ( 2014 ). We 
assume the AGB star to consist of a core with Mcore = 
0.535 Mq and an envelope with Mpn v = 1.465 M,^, con - 
sistent with the parameters inferred bv iMustill et al.l (120 13h . 
The mass enclosed within radius r is then given as 


M (r) — iWcore + 


Airr^ p{r)dr, 


(5) 


0.1 



1.0 10.0 
Radius 


100.0 


Fig. 1 Comparison of the gravitational binding energy of 
the envelope with the gravitational energy released via the 
inspiral in NN Ser. Both energies become equal for Eb = 
Eg ^ 1.2 X 10^^ erg at a r adius of 0.94 leading to 


the ejection of the envelope (ISchleicher & Dreizlea 1201 4ft . 


with Rmre- ^ 0-01 Rip - As in the model s of 

Nordhaus & Blackman ( 200^ . Tauris & DewiI ( 2001 ) and 


Soken (Il992h . we assume a power-law density profile in the 
envelope of the AGB star where 


p{r) = Ar 


(6) 


with w = 2. From the mass of the envelope, the normaliza¬ 
tion follows as 


A = 


M„ 


AnR^ 


(7) 


with i?* the radius of the AGB star, for which we adopt a 
typical value ii* ^ 185 R©. We consider now the inspi¬ 
ral of the secondary star with mass M 2 — 0.111 Mq and 
assume that the ejection occurs once that the released grav¬ 
itational energy exceeds the binding energy of the envelope 
mass outside radius r. The latter is given as 

rR, 


Eb = 


G{M{r) + M 2 ) 


p{r)dr, 


while the former is calculated via 

GM{r)M 2 GiM^ore + Me„v)M2 


Eg = 


2r 


2 i?* 


( 8 ) 


(9) 


In this equation, the factor 1/2 reflects that half of the gravi¬ 
tational potential is balanced by the kinetic energy. We com¬ 
pare both expressions in Fig.[T]and show that they are equal 


In the following, we will assume that the energy in¬ 
jection occurs almost instantaneously, which may be jus¬ 
tified due to the short timescales once that the secondary 
has reached small scales close to the central core. This is 
also the stage where most of the gravitational energy will 
be released. We can therefore employ the self-similar S edov 
solution for a power law density profile ( Sedovl 1959h . As 
the injected energy Eq, we consider the sum of the binding 
energy of the envelope Eb and the available thermal energy 
in the envelope, which we calculate from the virial theorem 
as 0.5Eb- The position of the shock front is then expected 
to evolve as 

/ fp x2 \ ^/ 

= (^) ’ 

where the constant a can be calculated from energy conser¬ 
vation. It is then possible to obtain the post- shock profiles of 
the radial velocity, density and pres sure as (IKashi & Sokei , 


201 ItlSchleicher & Dreizlen l2014l) 


1 r 

2G 


P = 


P = 


37 — 1 f 

2l^7±ii;^8/(7-ri) ^ 4Ar-^X^ 
r ‘^(7 — 1 ) 

A 2(7 4 - 1) ^ 1 

3 


j.ui-2f2 (^ 3 ^ _ l'^2 


2 ^3 


( 11 ) 


( 12 ) 


(13) 
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where we assumed ^ = hj^ after the second equality and 
introduced the self-similar variable 


A = 



l/(5-:v) 


r/Rs- 


(14) 


From the consideration that F^therm + ^'kin = Eq, one can 
show that a = TT. We evaluate the velocity profile at the time 
where Rs = R*, and compare with the escape velocity of 
the system. 


Vesc{r) 


^ 2G'(Mcore + ^ 2 ) ^ ^ 


(15) 


As a result, we find that the mass on scales above 106 R© 
becomes unbound. The ejected mass which remains bound 
to the system corresponds to Mbound ^ 0.133 Mq or 140 
Jupiter masses. In reality, the bound fraction may even be 
higher if the ejection process is highly inhomogeneous. 

A relevant question concerns the angular momentum of 
the ejected material. The total angular momentum that is 
deposited in the envelope by the secondary can be estimated 
as 


.^dep = M2RS, 


GMi 


i?* 


(16) 


yielding about 1.2 x 10®^ erg cm^ s“^. It is however likely 
that the envelope was previously rotating, as angular mo¬ 
mentum can be exchanged with the secondary via tidal 
torques before the onset of the common envelope phas e 
dHurlev et all 120021: IhuI Il98ll: ISokeil 1 19951: IZahni l2008h . 


As shown by Bear & Sokei I 2 OIC ), the rotational velocity 
can reach /^ot = 45% of the breakup velocity. The angular 
momentum due to rotation of the envelope is thus given as 




Trot — 


/rot • 4:Trr^drp{r)r 




(17) 


yielding a contribution of Trot 2.6x 10®^ ergcm^ s“^.To 
determine the spatial extent of the disk, the main quantity of 
interest is however the specific angular momentum, which 
we parametrize as 


Tdisk 

ATdisk 


L. 


= aL 


dep 


M, 


(18) 


Both due to the inhomogeneous injection of the angu¬ 
lar momentum of the secondary, as well as due to the 
initial rotation of the envelope, one may expect that a 
value of a T, ^ 10 can be achieved (see discussion in 


Schleicher & Dreizleri 120141) . Similar enhancements of the 
specific angular momentum have been obs erved for instance 


in the gaseo us disks in the Red Rectangle (iBuiarrabal et al. 


20051120031). We also note that an independent analysis of 
the angular moment um constraints has been pursued by 
Bear & Soken (l2014ll . finding no relevant constraints for 
planet formation in NN Ser. 


3.2 Planet formation via gravitational instabilities 


In the following, we assume that the gas of the gravitation¬ 
ally bound material settles into a disk described by a power- 
law profile of the gas surface density, 

E(r) = So(^)”, (19) 


where rout denotes the outer radius of the disk, Eq the disk 
surface density at the outer radius, and n the power-law in¬ 
dex. We will in particular c onsider the case of a Mestel disk 
with n = 1 (lMestelLll963h . We note here that Eg and Tout 
can be calculated from the disk mass and angu lar momen¬ 
tum as outlined in Schleicher & Dreizlei ( 20141) . and the re¬ 
sulting surface density profiles are given in Fig.|2] 



Radius [AU] 


Fig. 2 The gas surfac e density as a function of rad ius for 
the models outlined by ISchleicher & Dreizlen (120141) . with 
power-law indices from n = 0 to n = 1.5.The models 
are norma lized to reproduce the obser ved planetary mass 
at 5.4 AU ( Schleicher & Dreizlei . 2014 ). 


We expect that the disk efficiently cools through the 
emission of the dust grains, and the Toomre Q parameter 
reaches as marginally stable state with 


Q 


ttGE 


1 . 


We therefore have 


( 20 ) 


Cs 


ttGE 

n 


( 21 ) 


We note that U is calculated assuming Kepler rotation with 
a central mass corresponding to the mass of the pres ent sys¬ 
tem. The disk heigh t follows as h{r) = Cs{r)/fl{r) ((Levin . 
2007HLodatoll20()7l) . 


The mass of the initial clumps fo rming by th e grav¬ 
itational instability is then given as ( Boleye^,^ 2010 t : 


___ lev 

Meru & BateLl2010H201 ll: [Rogers & Wadslevll2012h 


Mol = E(r)/j2(r). 


( 22 ) 


© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 


www.an-journal.org 













































































Astron. Nachr. / AN (2006) 


793 


These clumps may grow on the orbital timescale until they 
reach the gap openin g mass where the whole d isk is evapo¬ 
rated near their orbit (ILin & PapaloizouiIl986l) : 


Mf = Mci 




1/2 


(D 


1/2 


(23) 


V 0.3 /- 

Here acrit denotes the a-parameter for viscous dissipation 
for which fragmenta tion occurs. We adopt here a generic 
value of acrit ^ 0.3 (iGammi el 120011) . The resulting clump 
masses are given in Fig. [3 We note in particular that it is 
straightforward to reproduce the characteristic mass scale 
of the planets, while the expected position has a stronger 
dependence on the specific angular momentum. In addition, 
the potential effect o f migration has to be considered (e.g. 
Baruteau et al.L 201 lb . As a result of such migration, t he ob- 


served resonances of the planets (IBeuermann et al.L l2013l) 
may form on timescales of a few orbital periods. 



Fig. 3 The f inal clump masses as a functi on of radius for 
the models by Schleicher & Dreizleil ( 2014 ) in the absence 
of radiati ve feedback, with power-law slopes from n = 0 to 
n = 1.5 ( Schleicher & Dreizlei , 20141) . 


3.3 Radiative feedback 


Radiation from the central star can substantially heat up the 
disk and stabilize it against gravitatio nal instabilities in the 
interior. We adopt here the approach of IChiang & Goldreich 
(Il997h to estimate the temperature in the midplane of the 
disk as 




1/4 


(t) 


1/2 




(24) 


where T* is the effective temperature of the central star and 
r* its radius. We consider here the stage very early after 
the ejection, and therefore estimate the radius of the star to 
be comparable to the remain ing radius of the core (see e.g. 
Schleicher & Dreizleiil2014l) . The grazing angle 0 describes 
the angle at which light from the star strikes the disk, and 
is estimated as 0 0.4r*/r. For the temperature, we then 



Radius [AU] 


Fig. 4 The gas temperature in the presence of photoheat¬ 
ing from the star in the c ase of a Mestel disk, correspon ding 
to model C described by Schleicher & Dreizleil ( 20141) . We 
explore the effect of different effective temperature of the 
central star T*. The minimum temperature in the model is 
based on the assumption that the dis k should be marginally 
stable ( Schleicher & Dreizlei , 2014 ). 


adopt the minimum of Eq. |24]and the temperature expected 
in a marginally stable state {Q ~ 1). 

The resulting thermal evolution is given in Fig.|4] show¬ 
ing that radiative feedback can stabilize the disk on scales 
around a few AU. The latter may naturally explain the sta¬ 
bility of the disk in the interior, and could in fact provide 
a natural point to stop the migration, as it becomes easier 
to create a gap in the disk once it is gravitationally stable 
( Baruteau et al. . 2011). 


3.4 Comparison with the observed population 

While NN Ser is currently the system with the most de¬ 
tailed data and the highest-quality fits to the planetary 
orbits, Zorotovic & Schreibed ( 2013 ) have suggested 12 
PCEB-systems which may harbor such massive planets, 
some of which even indicate the presence of two plan¬ 
ets as in NN Ser. These systems include five detached 
systems with a hot subdwarf B (sdB), four detached sys¬ 
tems with a white dwarf (WD) and three cataclysmic vari¬ 
ables (CVs). Eor the system NSVS14256825, the two- 
planet solution was sho wn to be dynamically unstable 


stable dBeuermann et al 


(IWittenm ver et al.L 12013) , whil e the one-planet solution is 


201^. In thi s case, we therefore 


use the data of Beuermann et al. ( 2012h. while we ad opt the 
parameters inferred bv IZorotovic & Schreibeii(l2013h in the 
other cases. 

The comparison of the theoretical predictions with the 
observationally inferred masses is given in Eig. |5] We find 
that there is a population showing good agreement with the 
theoretical predictions, as well as an additional population 
with even higher masses. The latter may hint either at a dif¬ 
ferent origin of these planets, or the potential presence of 
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additional effects related to magnetic activity, which cannot 
be ruled out in some cases. 


4 Discussion and conclusions 


In this article, we have summarized the main arguments 
against a hrst-generation scenario for the planets in NN 
Ser, and presented a theoretical model explaining the forma¬ 
tion of planets from the ejecta during the common envelope 
phase. The model naturally explains the observed planetary 
masses in NN Ser, and is in good agreement with a num¬ 
ber of additional systems. In addition, there seems to be a 
population where the observationally inferred masses sig¬ 
nificantly exceed the theoretical predictions. It is interesting 
to note that in exactly these cases, it is difficult to justify a 
second- generation scenario d ue to angular momentum con¬ 
straints (IBear & SokenuM^ . The latter can potentially hint 
at a different planetary origin. We note that beyond a purely 
hrst- or second-generation origin, also hybrid scenarios are 
conceivable where the existing planets accrete additional 
mass from the bound material remaining after the ejection. 



0 2 4 6 8 

Mi,obs / 


Fig. 5 Comparison of model predictons (diamonds) 
with planetary masses infe rred from observations by 
Zorotovic & Schreibeij ( 2()I^ (triangles). The data indicate 


the potential presence of two populations, one comparable 
to our model predictions, and one with even higher masses 
( Schleicher & Dreizler . 2014 ). 


The scenario proposed here can be validated through 
the further observation of the eclipsing binary systems, as a 
strict periodicity is expected in the case of planetary orbits. 
In addition, complementary approaches like the detection of 
the planets via their thermal emission could provide an im¬ 
portant pathway to confirm their existence. Telescopes like 
Gemini and ALMA may in addition look for remants of the 
gravitationally bound gas in NN Ser to identify a potential 
debris disk or the gas ejected during the common envelope 
event. The latter may provide further information to signih- 
cantly rehne such a model. 
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